Are Bosonic Replicas Faulty? 
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Motivated by the ongoing discussion about a seeming asymmetry in the performance of fermionic 
and bosonic replicas, we present an exact, nonperturbative approach to both fermionic and bosonic 
zero-dimensional replica field theories belonging to the broadly interpreted /3 = 2 Dyson symmetry 
class. We then utilise the formalism developed to demonstrate that the bosonic replicas do correctly 
reproduce the microscopic spectral density in the QCD-inspired chiral Gaussian unitary ensemble. 
This disproves the myth that the bosonic replica field theories are intrinsically faulty. 
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Introduction. — Since the mid-1990s, there has been a 
revived interest in the field theoretic approaches tailor- 
made to the analysis of interacting disordered and quan- 
tum chaotic systems. In particular, the exact Keldysh 
[H and approximate supersymmetry Q techniques have 
been conceived to offer a nonperturbative alternative to 
the notoriously known replica field theories [H, 13, HI whose 
legitimacy has been questioned [6[ for more than two 
decades. Sadly, the newly proposed field theoretic ap- 
proaches [l|, [3] have not yet evolved into efficient calcula- 
tional tools, and their success has been very limited. 

At the same time, substantial progress [1, 0, was 
achieved over the past few years in resolving controversies 
surrounding nonlinear replica a models. Specifically, the 
fermionic version [J, |ll| of a replica field theory consid- 
ered in the so-called zero-dimensional (random-matrix- 
theory [lij]) limit was proven [§| to be exactly integrable. 
This observation brought into play the whole machinery 
of the theory of integrable hierarchies [3] and eventu- 
ally resulted in reconstructing [U, Q the exact spectral 
densities and/or correlation functions for the paradig- 
matic Gaussian unitary ensemble (GUE), the QCD in- 
spired chiral GUE (chGUE) [3] and Ginibre's ensemble 
[151 ] of complex non-Hermitean random matrices. The ex- 
act fermionic replicas also have immediate implications 
for the nonperturbative physics of the ID lattice impen- 
etrable bosons 16]. The supersymmetric variation [101 ] 
of exact replicas [8( has already produced important new 
results 17| for the QCD at nonzero chemical potential. 

The present Letter, prompted by the ongoing dis- 
cussion about a seeming asymmetry in 



the performance of fermionic and bosonic replicas [21 1 
(fermionic-bosonic dichotomy) , addresses the problem of 
integrability of zero-dimensional bosonic replica field the- 
ories much in line with the ideas of Refs. [1, Hav- 
ing formulated a general nonperturbative theory of both 
fermionic and bosonic replicas, we further concentrate 
on the chGUE matrix model and develop an integrable 
theory of the corresponding nonlinear bosonic replica a 
model. Contrary to the claims made in the literature 
[la ], the latter is shown to produce the exact expres- 
sion for the chGUE density of eigenlevels in the physi- 



cally relevant limit of infinite-dimensional matrices. This 
achievement, representing the main outcome of our study, 
provides strong evidence that the bosonic replicas are as 
good and reliable as the fermionic ones. We conjecture 
that the above statement holds in the whole generality, 
no matter what particular random matrix model is being 
treated. 

How replicas arise and why they are tricky. — Replica 
field theories (be it the original bosonic formulation in- 
vented by Wegner Q or its fermionic counterpart Q fur- 
ther extended by Finkelstein Q to accommodate the in- 
teraction effects) are based on the identity 



log Z = lim 



Z n - 1 



(1) 



which can be very useful 22j in evaluating the average 
(logZ). Upon assigning Z the meaning of a quantum 
partition function Z{<;) = YY a -i det(? a — "K) of a system 
characterised by a stochastic Hamiltonian JC, the identity 
(P) can be utilised to represent the average p-point Green 
function G(«j) = (IIa=i Tr (<«* — IK) 1 ) in terms of the 
average characteristic polynomials 

Z n (s) = (j[<let n (<; a -X)\, ? = (&,••• ,Sp), (2) 

to be referred to as the replica partition function. Notice 
that Eq. © is defined for n S R 23]. The recipe, known 
as the replica limit, reads 

G(s)= lim ±-^...3 Zn {s). (3) 

Equation ([3]) assumes a mutual commutativity of the fol- 
lowing operations: the replica limit, differentiation, dis- 
order averaging denoted by the angular brackets (■••), 
and a thermodynamic limit, if necessary. 

Seemingly innocent at first glance, the prescription (J3]) 
is much trickier than one could naively expect. Indeed, 
in order to calculate the replica partition function Z n (?) 
nonperturbatively, a field theorist interprets Z n (<;) in 
Eq. f2]) as a substitute for |n| € Z + identical noninter- 
acting copies, or replicas, of the original random system. 
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Each copy, exemplified by the product n«=i det^Q, — Di) 
of p > 1 single determinants, is represented by a func- 
tional integral over an auxiliary field which is either 
fermionic or bosonic by nature, depending on the sign 
of n. Exponentiating a random Hamiltonian !H, such 
a representation facilitates a nonperturbative averaging 
over the ensemble of stochastic Hamiltonians in Eq. ([2]) 
and eventually results in effective field theories defined 
on either a compact Q (fermionic, n £ Z + ) or a non- 
compact [H (bosonic, n £ Z~) manifold. Such a replica 



mapping, Z„ eR (<?) 



Z„ g z±(Oj clearly indicates the 



key problem of replicas. By derivation, the validity of 
Z„ gZ ±(^) is restricted to n £ Z^, which is not enough 
for implementing the replica limit (|3|) determined by the 
behaviour of Z n (s) in the vicinity of n — 0. This mis- 
match between the "available" and the "needed" is at the 
heart of the trickery with which the replica field theories 
are often charged [6J]. 

The canonical way to bridge this gap is to determine 
Z n (<;) for n £ TL^ , and then attempt to analytically con- 
tinue Z n ($) away from n integers, in general, and to a 
proper vicinity of n = 0, in particular. Since perform- 
ing an analytic continuation based on an approximate 
result is a mathematically questionable procedure, the 
evaluation of Z n (s) must be done exactly. Below, we 
will show how such a nonperturbative calculation can be 
carried out in quite a general setting. The approach to 
be presented applies to the matrix models belonging to 
the broadly interpreted Dyson's (3 = 2 symmetry class 
[H, HH and is by far more flexible and efficient than the 
oneofRef. Q. 

Nonperturbative approach to replicas. — Let us concen- 
trate on the fermionic and/or bosonic replica field theo- 
ries whose mapped partition functions admit the eigen- 
value representation (n is supposed to be positive) 

« n 

^ f/b) W=/ Y[d\ k T(<;;\ k )e- v ^^Al(\). (4) 



fc=l 



Here V n (A) is a "confinement potential" which may de- 
pend on the replica index ±n; T(s; A) is a function accom- 
modating relevant physical parameters <r of the theory 



[they are not necessarily the energies specified in Eq. ©J. 
To treat the fermionic and bosonic replicas on the same 
footing, the integration domain T> was chosen to be [25| 
23 = [fj=i[c2j-i,c 2j ]. 

To detemine the replica partition function Zn^ h \<i) 
nonperturbatively, we adopt the "deform-and-study" ap- 
proach, a standard string theory method of revealing hid- 
den structures. Its main idea consists of "embedding" 
^ f/b) (^) into a more general theory of r functions 



r(») 



i r - 



x e 



(5) 



which posses the infinite-dimensional parameter space 
t = (£j.,t2( ' • ' ) arising as the result of the ^-deformation 
v(t; A) = ^j—itjX ; the auxiliary parameter s is as- 
sumed to be an integer, s £ Z. Studying the evolution 
of t functions in the extended (n,s,t, space allows 
us to identify the highly nontrivial, nonlinear differen- 
tial hierarchical relations between them. Miraculously, a 
projection of these relations, taken at s = 0, onto the 
hyperplane t = 0, 



t = 

«=0 



(6) 



generates, among others, a closed nonlinear differen- 
tial equation for the replica partition function Zn b \$). 
Since this nonperturbative equation appears to contain 
the replica (or hierarchy) index n as a parameter, it is 
expected [8[ to serve as a proper starting point for build- 
ing a consistent analytic continuation of Zn^ h \<s) away 
from n integers. 

Having formulated the crux of the method, let us turn 
to its detailed exposition. The two key ingredients of the 
exact theory of r functions are (i) the bilinear identity 
[IH and (ii) the (linear) Virasoro constraints [26] . 

(i) The bilinear identity encodes an infinite set of hi- 
erarchically structured nonlinear differential equations in 
the variables {tj}. For the model introduced in Eq. 
the bilinear identity reads 27, 28|: 



, HI 



Here, a £ R is a free parameter; the integration contour 
Coo encompasses the point z = oo; the notation t± [z^ 1 ] 
stands for the infinite set of parameters {tj ± z~ J /j}; for 
brevity, the physical parameters <r were dropped from the 
arguments of r functions. 

Being expanded in terms of t' — t and a, Eq. (JT]) 



generates four integrable hierarchies. One of them, the 
Kadomtsev-Petviashvili (KP) hierarchy in the Hirota 
form [2j| 

i^ fc r«(t)o T W(t) = afc+1 ([r>])r«(t)o T W(t) (8) 
(k > 3) is of primary importance for the exact theory 
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of replicas 3l|. The first nontrivial member of the KP 
hierarchy reads 



+ 6 (d 2 tl logrWfct)) 2 = 0. 



(9) 



In what follows, it will be shown that its projection onto 
s = and t = [Eq. J^J] gives rise to a nonlinear differen- 
tial equation for the replica partition function Zn^ h \<;). 
(ii) Since we are interested in deriving a differential 

equation for Zn (?) in terms of the derivatives over 
physical parameters we have to seek an additional 
block of the theory that would make a link between the 
{tj} derivatives in Eq. ([9]) taken at t = and the deriva- 
tives over physical parameters {?j}. The study [321 ] by 
Adler, Shiota, and van Moerbeke suggests that the miss- 
ing block is the Virasoro constraints which reflect the in- 
variance of the r function [Eq. ©J under a change of the 
integration variables. In the present context, it is useful 
to demand the invariance under the transformation 

dim(c ) 

Aj-Mj+eMfVfoi) II C/*j-<£). e>0 < ( 10 ) 
fc=i 



where c' = {c%, ■ ■ ■ , C2 r } \ {±00} ■ The function /(A) is 
related to the confinement potential V n - s (X) through the 
parameterisation dV n - s /d\ = —g(X)/f(X), where g(X) = 
YlJLo bjX^ and /(A) = Y^jLa % A J depend on n and s. 

The transformation (fT0|) induces the Virasoro-like con- 
straints that can be written in the additive form 



rW(?;t) = 0, 4 > -1 



(11) 



The operators H%(t) and £jq(t) are associated with the 
e -v„_ s (A) an( ^ ^. ne r(^; A) parts of the integrand, respec- 
tively; also, £g =1 (t) = 0. The first operator &Y(t) can 
be expressed in terms of the Virasoro operators [26| 



L q {t) 



3*3 9ta 



9 

E 

J=0 



dud t . 



?>-l, (12) 



which depend solely on the symmetry of the replica field 
theory and obey, for all p, q > —1, the Virasoro algebra 
[£jp,L q ] = (p— q)L p+q . [Equation (fT2|) assumes that dt 
is identified with the multiplicity of the matrix integral 
in Eq. ©, d to = n]. Explicitly, it holds ^ that 



dim(c') 00 

&q(t) = ^ s dim(c')-fc(-[ <T ])E ( a e&q+k+i(t) - b t d t 

£=0 



k=0 



q + k + t+l I ) 



dim(c') 

Ai = ~. E W' ?>-l. 



(13) 



k=l 



Here, Sfe(i) are the Schur polynomials [301 ] . 

While very similar in spirit, the calculation of L^{t), 
the second ingredient in Eq. (jlip . is more of an art since 
the function A) in Eq. ([5]) may significantly vary from 
one replica model to the other. 

Remarkably, for t — 0, the two equations [Eqs. ([9]) 
and (fTT|) ] can be solved jointly to bring a closed non- 
linear differential equation for Z„ h \<;). It is this equa- 
tion which, being supplemented by appropriate boundary 
conditions, provides a truly nonperturbative description 
of the replica partition functions and facilitates perform- 
ing the replica limit. 

Chiral GUE and bosonic replicas. — To see the above 
formalism at work and also answer the question raised in 
the title of our Letter, let us consider the N x N chGUE 
random matrices 



W 

wt 



(14) 



known to describe the low-energy sector of SU(7V C > 3) 
QCD in the fundamental representation [14|. Composed 
of rectangular ul X ur random matrices W with the 



Gaussian distributed complex- valued entries 



2tt 



exp 



iVE 2 



trW f W 



,(15) 



where N = til + tir, the matrix !Kd has exactly 
v = \ur — nj_\ zero eigenvalues identified with the topo- 
logical charge v] the remaining eigenvalues occur in pairs 
{±Aj}; the parameter £ denotes the chiral condensate. 

To determine the (microscopic) spectral density from 
the bosonic replicas, we define the replica partition func- 
tion zji\<;) = (det~"(<; + i!Kd)) w and map it onto a 
bosonic field theory. In the half-plane Re<; > 0, the par- 
tition function Z^ (<;) reduces to [3, [2(| 



dn n (Q) det^- n O exp — ^Tr(Q + Q -1 ) 



(16) 

where the integration domain §„ spans all n x n posi- 
tive definite Hermitean matrices Q. Equation flU]) was 
derived in the thermodynamic limit N — > 00 with the 
spectral parameter oj = ^iVE being kept fixed (Reuj > 0). 

Spotting the invariance of the integrand in Eq. (fT6|) 
under the unitary rotation of the matrix Q, one readily 
realises that Z^\uj) belongs to the class of r functions 
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specified by Eq. ([5]) where D is set to R + , the potential placed with T{ui; A) = exp [— (w/2)(A + A 1 )] . This ob- 
V n -s is V n - S {\) = {n-s-v) log A, and rfc; A) is re- ser vation implies that the associated r function T^ s \uj;t) 

I satisfies both the first KP equation ([9]) and the Virasoro 

constraints JTTJ) with [53] 

iV(t)=i q+1 (t) + (v-n + s)d tq+1 , t r q (u ! ;t) = -^d tq+2 -6 g ,- 1 (u ) d u , + ^d tl )+[l-5 q! - 1 } ^d tq . (17) 
Projecting Eq. ([9]) taken at s = onto t = 0, and expressing the partial derivatives therein via the derivatives over u 



with the help of Eqs. (fTTj) and (fTT|) , we conclude that 2T„ (w) — n\ Tn\to; 0) obeys the differential equation [27 1 



K + - (4 + 1 + 4( 1 + ^ 2) ) + 6(K) 2 + 1 - 4( " 2 3 + - 2) ^ - + A -K K + A 4 = o (18) 



that can be reduced to the Painleve III. Here 
h n {oj) — duj \ogZ^n\uj). Considered together with 
the boundary conditions h n (u> — > 0) ~ —nvjio and 
/i n (w — > 00) ~ — n — n 2 /(2u), following from Eq. |5]), 
the nonlinear differential equation (|18p provides a non- 
perturbative characterisation of Z„ (lj) for all n 6 Z + . 

To pave the way for the replica calculation of the 
Green function G(u>) determined by the replica limit 
G(uj) — — lim„^o n ~ 1 h n (uj) 1 one has to analytically 
continue hJUS) away from n integers. The previous 
studies d, Q suggest that the sought analytic contin- 
uation is given by the very same Eq. (fl8|) where the 
replica parameter n is let to explore the entire real axis. 
This leap makes the rest of the calculation straightfor- 
ward. Representing h n (u>) in the vicinity of n = as 
h n (u>) — YlpLi nFa p{ UJ )i we conclude that G(u>) — 
— <Xi(u>) satisfies the equation 

uj 3 G"' + 2u?G" - (l + 4v 2 + Aoj 2 ) ujG' + (1 - Av 2 ) G = 0. 

(19) 

Its solution, subject to the boundary conditions consis- 
tent with those specified below Eq. ([T5]). brings the mi- 
croscopic spectral density q(u) = 7r _1 ReG(iw-|-0) in the 
form 



q(u) = vS(uj) + ^ Jl{u) - J v -\(<J)Jv+i(u) 



(20) 



Obtained within the framework of bosonic replicas, this 
celebrated formula provides strong evidence against the 
idea of their inapplicability to the nonperturbative de- 
scription of random matrix spectra [fj, in general, and 
of the chGUE spectra [l8( , in particular. In view of the 
previous study [8j on the performance of fermionic repli- 
cas, we are led to speculate that truly nonperturbative 
approaches to nonlinear replica a models leave no room 
for the fermionic-bosonic dichotomy. 
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